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Abstract Metallic behavior of band-insulator/ Mott-insulator interfaces was ob- 
served in artificial perovskite superlattices such as in nanoscale SrTiOs/LaTiOa 
multilayers. Applying a semiclassical perspective to the parallel electronic transport 
we identify two major ingredients relevant for such systems: i) the quantum con- 
finement of the conduction electrons (superlattice modulation) leads to a complex, 
quasi-two dimensional subband structure with both hole- and electron-like Fermi 
surfaces, ii) strong electron-electron interaction requires a substantial renormaliza- 
tion of the quasi-particle dispersion. We characterize this renormalization by two 
sets of parameters, namely, the quasi-particle weight and the induced particle-hole 
asymmetry of each partially filled subband. In our study, the quasi-particle disper- 
sion is calculated self-consistently as function of microscopic parameters using the 
slave-boson mean-field approximation introduced by Kotliar and Ruckenstein. We 
discuss the consequences of strong local correlations on the normal-state free-carrier 
response in the optical conductivity and on the thermoelectric effects. 



1 Introduction 

Recent experiments O] have shown that a metallic state can be stabilized at the 
interface between the Mott insulator LaTiOs and the band insulator SrTiOs. There 
is strong evidence that in such systems electronic charge is redistributed between 
the Mott insulator (MI) and the band insulator (BI) in order to compensate for the 
mismatch of the work functions and to avoid the so-called polar catastrophe (|2l. The 
electronic charge reconstruction Q at the interface leads to metallic behavior. 

The experimental data H] |4] |5] |6l are consistent with Fermi liquid behavior and 
a single-particle perspective, where transport properties are studied by the semiclas- 
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sical transport equations, offers a natural starting point. However, it is necessary 
to clarify how the single-particle picture for weakly interacting electrons, which 
is successfully applied in the study of semiconductor nano-structures, is modified 
by strong electronic correlations. In particular, understanding the renormalization 
of the quasi-particle dispersion as function of microscopic parameters is crucial. 
Let us in the following ignore complicating aspects related to the orbital degrees 
of freedom p] or to possible symmetry -broken phases lO. Then, from quite gen- 
eral considerations, we can expect two major ingredients determining the electronic 
structure; 

(i) The LaTiOs bulk system has a Ti-3c/' configuration whereas the SrTiOs com- 
pound has a Ti-3d^ configuration. Therefore, from a single-particle point of view, 
the BI/MI/BI sandwich acts as a quantum well confining the conduction electrons 
to the MI region, cf. Fig. [T| The bound states of the quantum well form quasi-two 
dimensional subbands with dispersion labeled by the in-plane momentum k 
and the subband index vQThe Fermi surface (FS) defined by the k points satisfying 
£kv = contains in general both open and closed sheets which brings about that 
electron-like and hole-like contributions can lead to partial compensation lITOl . 

(ii) Strong electron- electron interaction introduces novel electronic physics at 
the band-insulator/ Mott-insulator interface. When the local self-energy corrections 
are dominant, we can assume that istv = £'v(ek), where is the non-interacting 
in-plane dispersion |l8]|9l - In this case, the renormalization of the quasi-particle dis- 
persion is characterized by two sets of parameters. On the one hand, the on-site 
repulsion leads to a reduction of the Fermi velocity of the subband v by a factor 



(1) 

FS 



which is equal to the quasi-particle weight of the subband v. On the other hand, 
at the interface, the hybridization of (almost) localized with itinerant degrees of 
freedom induces an enhanced particle-hole asymmetry. We quantify this asymmetry 
by the dimensionless parameter 



'\ I ( '^^^^ 
del ') V ^e^. 



(2) 

FS 



where we have defined ZiCk = — fif, with Ej, the energy of the lower band edge. 
Starting from a microscopic model, we calculate the quasi-particle dispersion in a 
self-consistent way to obtain the dependence on microscopic parameters. We obtain 
an interfacial heavy-fermion state and discuss correlation effects as characterized 
by Eq. ([T]l and cf. Fig.|2] Implications for transport are illustrated by calculating 
the free carrier response and the thermoelectric effects, cf. Fig. [3] We find that in 
both cases important contributions arise from the interface. 



Our discussion is restricted to large superlattice periods where the dependence on the perpendic- 
ular momentum Q can be neglected when studying the parallel transport. 
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Fig. 1 The charge density 
Hi, the fraction of doubly 
occupied sites df and the 
Lagrange multiplier A; as 
obtained by the present mean- 
field approach for a band- 
insulator (BI)/Mott-insulator 
(MI) superlattice. L = N + 
M denotes the number of 
unit cells of the superlattice 
modulation where N is the 
number of Ml-layers and M 
the number of Bl-layers. The 
set of microscopic parameters 
is given in the figure. 



N = M = 10, Ur = 22t, Ec = O.St 
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2 Microscopic model 

For the microscopic description we assume perfect lattice match between the two 
materials, thereby neglecting aspects related to the lattice relaxation ifTTI . The mi- 
croscopic model is given by an extended single-orbital Hubbard model on a cubic 
lattice (introduced in Ref. Tl3l ) 

H = H,+Hu+Hee+ Hei + Hii. (3) 

Here, the kinetic energy is given by a nearest-neighbor tight-binding model and the 
on-site repulsion is modeled by a Hubbard interaction, 

H, = -t c\^Cj„ ^h.c. and Hu ^ Uj^rii^tiii, (4) 

where iiia = cJ^-Cj^. The nanoscale structure is defined by the superlattice period L = 
N + M and the number of counter-ion layers (see also Fig.[T]). They simulate the 
difference between Sr"+ and La^+ and sit in the center between the electronic sites 
||3l interacting with the electrons through the long-range electron-ion interaction 

Hei = -Ec y - — (5) 

where r'°" denotes the position of the ions, «, = «,| + n,j, and we have introduced 
the parameter Ec controlling the screening length. Furthermore, the long-range 
electron-electron and ion-ion interaction energies are given by 
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respectively. The number of electrons is fixed by the charge-neutraUty condition. 
Notice that we can formally relate the parameter Ec to an effective dielectric con- 
stant e = /Eca where a is the lattice constant and e > the elementary charge. 
However, for a more realistic description of the screening at the interface a single 
parameter for the long-range electron-electron interaction is too crude. In fact, the 
polarization of the lattice dominates the dielectric constant in the considered tran- 
sition metal oxides and the effect of the relaxation of the lattice near the interface 
introduces additional parameters in an effective model description lflTl[T2l . For sim- 
plicity, such effects are not considered here. 



3 Slave-boson mean field approximation 

To discuss the low-energy behavior of the model (|3]l in the normal state we ap- 
ply the four-boson mean-field approximation of Kotliar and Ruckenstein il4|. This 
approach allows to discuss the effect of local self-energy corrections by introduc- 
ing auxiliary bosons representing the local charge and spin degrees of freedom to- 
gether with pseudo fermions. The effective low-energy theory is then obtained by 
the saddle-point approximation for the slave bosons which can be controlled by 
a 1 /A^-expansion of a suitable generalization of the slave-boson action ifTSl [161 . 
The remaining fermionic degrees of freedom are interpreted as the Landau quasi- 
particles of a Fermi liquid which are dressed by the interactions and therefore have 
modified single-particle properties. 



3.1 Superlattice geometry 

For a quantum well system, the effective low-energy model was derived in Ref. HI. 
In the following, we will briefly discuss the case of a superlattice. Assuming a trans- 
lational invariant state in the in-plane direction, the problem of finding the eigenval- 
ues of the effective low-energy Hamiltonian reduces to a one-dimensional problem 
which is parameterized by the non-interacting dispersion Ek = — 2f (cos^^^ + cosfc,.); 



Here, I labels the layers along the direction of the superlattice modulation. The hop- 
ping renormalization amplitude zi depends on the charge density «/ and the fraction 
of doubly occupied sites df in layer / and is given by the standard Gutzwiller ex- 
pression ifm 



(z^ek + A/)v/(/)-f zizi+yW{l + y) = E\if{l)- 

y=±l 



(7) 




(8) 
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The mean fields «/ and di are determined by the minimum of an appropriate free 
energy, as discussed below. The Lagrange multiplier Xi acts as a single-particle po- 
tential and enforces the self-consistency of the electronic charge distribution. 

Decomposing Xj/ according to ^fKv{l)e'^', where K = {k,Q), ~n/La < 

Q < n/La with i/kv(/ +L) = y/Kv{l), the problem reduces to diagonalizing the 
following matrix 

/ zfek + Ai -tzizie'O ... -fziZLe^'^X 
-tZ2Zie-'Q zje^ + h -tZ2Z3e'Q ... 

V -tZLZie'Q ... -tZL-iZte-'^ zle^ + h J 

A further simplification is obtained by restricting to superlattices with a large pe- 
riod L 1 . In this case, the Q-dependence can be safely neglected for the parallel 
transport and it is sufficient to consider only K{Q). In this case, the quasi-particle 
dispersion has the form stated in the introduction, E ^Ey (et), V = 1 , . . . ,L. 



K{Q) = 



3.1.1 Long range Coulomb interaction 

The long-range Coulomb interaction is treated in the Hartree-type of mean-field 
calculation IS) . In order to find the interaction matrix W/// between electrons in layer 
/ and /' of the superlattice unit cell we explicitly take into account the periodicity of 
the charge distribution. In the same way, we also determine the resulting (screened) 
potential V/ of the counter-ions. In the end, these potentials have to be found self- 
consistently by simultaneously solving the mean-field equations and the Poisson 
equation, A(j){z) = —47t^j^, where <j){z) is the electrostatic potential and p(z) the 
charge distribution. 

We start with the approximation commonly found in the literature ifTSl . namely, 
we replace each layer by a uniformly charged plane, thereby respecting the polar 
nature of the Mott insulator. Corrections due to the discrete nature of the charge 
distribution (lattice) are calculated numerically, but are only significant very close to 
(or within) the considered layer lfT9l . It is convenient to use the following elementary 
solution (j)o determined by a periodical array of uniformly charged layers with period 
La 

Poiz) = L Siz + mLa)-p (10) 

mez 

where <7 = e/a^ is an elementary surface-charge density and p = a /La is a uniform 
background charge to keep the total system charge neutral. The solution can be 
written in a compact form by use of the polylogarithm Lii(z) = Ln>o 



en 



(11) 
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Between two neighboring layers at nLa and (n + 1 )La, the resulting potential is 
simply given by the parabola 

^o{z) = —-^{z — nLa){nLa+La~z)- (12) 

J-iC- 

From the elementary solution Eq. ( fTTT i it is straightforward to determine W/// and 
Vi by summing up the contributions from the different layers within the superlat- 
tice unit cell and adding the numerically determined correction terms due to the 
discreteness of the charge distribution. 



3.1.2 Free energy 

Eventually, we give the expression for the free energy per lattice site (j3^' = ksT) 

/(„,^,A) = -^^ln(l+.-/^^''v) + ^^rf,2^1^„,W,,«,-i^ 
P'^lr kv / zi. I 

(13) 

Sums over layers are restricted to a single superlattice unit cell. The self-consistency 
equations are solved by maximizing Eq. iT3[ with respect to A; and minimizing the 
resulting function with respect to the mean fields «/ and di under the constraint of 
charge neutrality, Y.I = 

A typical solution of the self-consistency equations at T = is shown in Fig. [T] 
for a superlattice with N = M = 10. The charge distribution «/ allows to naturally 
distinguish between three different regions, «/ « 1, «/ w 0.5 and «/ « 0. The lay- 
ers with filling «/ w 0.5 seperate the "Mott-insulating" (MI) regions («/ « 1) from 
the "band-insulating" (BI) regions with «/ « (notice that the whole system is ac- 
tually metallic). The presence of the interface (IF) layers with filhng n/ « 0.5 is a 
consequence of the polar nature of the Mott insulator. 



4 Transport properties 

From the self-consistent solution of the mean-field equations we also obtain the 
quasi-particle dispersion ^tv and the envelope wave-function V^kvO- At low tem- 
peratures, it is expected that the transport properties can be understood from the 
properties of the quasi-particles lH. We therefore start by discussing the electronic 
structure. 
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Fig. 2 a) The quasi-particle 
dispersion Sky as function of 
the non-interacting in-plane 
dispersion e^. b) The sub- 
band filling iiy, c) the quasi- 
particle weight Zy and d) the 
induced particle-hole asym- 
metry Uy for the individual 
subbands which correspond 
to the dispersion in a). Circles 
(green) are associated with the 
subbands of the MI region, 
squares (red) with the IF and 
diamonds (blue) with the BI 
region. 



N = 15, M = 5,Ur = 22t. Ec = 0.8i 

2r 
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4.1 Characterization of the generic electronic structure 

The quasi-particle dispersion is shown in Fig.|2]a) for the = 15, M = 5 super- 
lattice. The value of the on-site repulsion, U,- = 22t, is well above the critical inter- 
action strength of the Mott transition in the half-filled bulk system ||20| , Uc « 16f . 
Thus, the quasi-particle dispersion shown in Fig.|2]corresponds to the strongly cor- 
related regime. In panels b) to d) we show the subband filling ny, the quasi-particle 
weight Zy [Eq. ([T]i] and the induced particle-hole asymmetry ay [Eq. ©J of the 
subbands shown in a). Although the spatial weight of the envelope wave-function 
V/kv(0 extends over the whole super unit cell and also depends on the value of 
£k 111 J it is possible to group the different subband states according to the regions 
where most of their spatial weight is located. We thus define Vmi = 1, . . . — 1, 
ViF = N,N +1 and Vbi ^ N + 2, . . . ,L. Notice that due to strong local correlations, 
the quasi-particle weight Zy^, and Zyjp is strongly reduced for the subbands of the 
MI and IF region [panel c)], whereas the particle-hole asymmetry ayjp is enhanced 
most dominantly for the subbands of the IF region [panel d)]. These are the basic 
characteristics of the interfacial heavy-fermion state obtained by the slave-boson 
mean-field approximation at low temperatures. The coherent hybridization of the 
itinerant degrees of freedom in the BI and IF region with the almost localized de- 
grees of freedom in the MI region is mediated by the intra-layer hopping and leads 
to heavy-fermion behavior of the interfacial subbands. The situation is reminiscent 
of heavy-fermion systems as described for example by the periodic Anderson model 
II2TII . However, in the present case, localized and itinerant degrees of freedom have 
the same orbital character but are separated spatially - in contrast to the classical 
heavy-fermion systems where the localized f-electrons hybridize with the states of 
the conduction band. 
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Fig. 3 (color online), (a) The Drude weight D as function of N jL for different values of Ec for a 
superlattice with period L = 20. (b) The total Seebeck coefficient S (triangles, black) and the dif- 
ferent contributions S„ associated with the MI (circles, green), IF (squares, red) and BI (diamonds, 
blue) regions as a function of Ec evaluated for a constant relaxation time. In both panels, the value 
of the on-site interaction is fixed at Ur = 22f . 



4.2 Drude weight 

The Drude weight in the optical conductivity is obtained as described in ||8]|9]. One 
finds the famiUar expression for a quasi-particle response, 

D = ^D., Dv = 4^A.v., (14) 

where Ay is the Fermi surface volume of the sheet v and Vy = ZydV^EkDrs/^ is 
the Fermi velocity averaged over the Fermi surface. In Fig.|3]a) we illustrate how D 
evolves from the band insulator (A^ = 0) to the Mott insulator (A^ = L), which both 
have a vanishing D. The maximal D as a function of the averaged electronic density 
N /L depends on the value of Ec/t and shifts to lower A^'s for increasing Ec because 
the screening length is reduced. 



4.3 Seebeck coefficient 

Thermoelectric effects are characterized by the Seebeck coefficient S. In terms of 
the subband contributions it is written as 

5=1^, <y = Y.<yv (15) 

y y 

Here, a is the total electrical conductivity and Sv the Seebeck coefficient associated 
with the subband v. In lowest order in the temperature we have ID 
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3 e ZvAe^ 



(16) 



where zie* = e* — e^, is measured from the band edge and defined through ^vlev) ~ 
0, Tv(e) is the relaxation time of the subband v and 

-^v(e) = |^02l^kek|'5(e-ek). (17) 

The prime (') in Eq. ( fTSI l denotes the derivative with respect to e at the Fermi sur- 
face. Attempts to calculate Ty from a microscopic model offers a challenging task in 
correlated, disordered and inhomogeneous systems (see also Ref. ID). For simplic- 
ity, we discuss here the case of a constant relaxation time. We assume an energy- 
independent relaxation time T„ for the subbands associated with the different re- 
gions, n — MI, IF, BI, and therefore obtain 

5 = ^^, 5„ = ^^, A, = ^Z)v, a„=D„T„. (18) 

Figure[3]b) shows the different contributions 5„ to the total Seebeck coefficient for a 
superlattice with = 15 and M = 5 as a function of the parameter Ec ■ The contribu- 
tion Sip from the subbands associated with the interface is most dominant and \Sip\ 
increases for increasing Ec (remember that a large value of Ec yields a sharp charge 
distribution). This can be understood by the fact that the particle-hole asymmetry 
Uv induced by the correlations is largest for the subbands of the interface region and 
increases for a sharper interface due to a reduction in the hybridization (intra-layer 
hopping). However, the total (absolute) Seebeck coefficient is at best equal to its 
largest subband contribution, \S\ < maxy \Sv\, but is in general smaller, as shown in 
Fig.[3]b) by assuming Tmi = Tip — Tbi. Nevertheless, we can state the conditions for 
which the interface contribution becomes large in the constant-relaxation-time ap- 
proximation: (i) strong electronic correlations, Ur > Uc- (ii) large values of A^, such 
that bulk-like properties in the center of the MI are obtained, (iii) a sharp interface 
(Ec > t). 



4.4 Comparison to the atomic limit result 

Note that the above discussed enhancement of |5if| at low temperatures is a non- 
local effect. It is therefore expected that any reduction of spatial coherence across 
the interface will suppress and eventually destroy this mechanism. The influence of 
the reduction of spatial coherence due to thermal fluctuations is clearly seen by an- 
alyzing S in the high-temperature (atomic) limit. Let us in the following discuss the 
situation t <^ ksT ^ U such that doubly occupied sites are completely suppressed. 
In the homogeneous system at density n, the thermopower S is given by the entropic 
contribution alone 1221123)1 . This yields Heikes formula 
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c 1 

i = log 



2(1 -«) 



(19) 



where the log 2 contribution arises from the entropy of the spin degree of freedom. 
As pointed out in Ref. |[24l, in transition metal oxides also the inclusion of orbital 
degrees of freedom is necessary. This is in principle straight forward but we will not 
discuss it further here. For the itihomogeneous system an appropriate generaUzation 
of Eq. ( fT9] l is given by 



2(l-«f) 



ni{\-ni) 



Here we have used the local quantities 



Si = log 



2(1 -n/) 



^1 = — 



ni{\~ni), 



(20) 



(21) 



with a temperature and doping independent constant A 11251 . The weighted sum in 
Eq. ( |20] | clearly shows that an inhomogeneous system is not favorable as long as 
a local description is appropriate. In fact, if one seeks to optimize the powerfactor 
PF= S^a in this limit for a spatially varying density profile, the optimal solution is 
found to be the homogeneous solution with optimal density n k, 0. 12, c.f. Ref. 1251 . 
This is exactly the opposite behavior than found in the low-temperaure limit. Re- 
stricting to purely electronic contribution, we therefore conclude that a spatially 
non-uniform system can only be favorable if spatial coherence is sustained. 



5 Conclusions 

In summary, we have studied aspects of the parallel transport at low-temperatures in 
strongly-correlated superlattices from the semiclassical point of view. The generic 
electronic structure is discussed and self-consistently computed from microscopic 
parameters using the four-boson approach of Kotilar and Ruckenstein to deal with 
strong local correlations. Implications for the parallel transport are illustrated by the 
free-carrier response and the thermoelectric effects. The presence of the interface 
introduces new aspects not feasible in the bulk systems. Here, we have discussed 
the scenario of an interfacial heavy-fermion state where the coherent hybridization 
of itinerant and almost localized degrees of freedom leads to a large particle-hole 
asymmetry which can be responsible for a high value of the Seebeck coefficient. 
We find that this mechanism is a non-local effect and that spatial coherence is cru- 
cial. If such mechanisms are also relevant for the giant thermopower observed in 
SrTi03/SrTio.8Nbo.203 superlattices ||26ll need to be clarified by further studies of 
more realistic models. 
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